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Abstract 

Background: Progressive remodeling of tine left ventricle (LV) following myocardial infarction (Ml) can lead to 
congestive heart failure, but the underlying initiation factors remain poorly defined. The objective of this study, 
accordingly, was to determine the key factors and elucidate the regulatory mechanisms of LV remodeling using 
integrated computational and experimental approaches. 

Results: By examining the extracellular matrix (ECM) gene expression and plasma analyte levels in C57/BL6J mice 
LV post-MI and ECM gene responses to transforming growth factor (TGF-pi) in cultured cardiac fibroblasts, we 
found that key factors in LV remodeling included macrophages, fibroblasts, transforming growth factor-pi, matrix 
metalloproteinase-9 (MMP-9), and specific collagen subtypes. We established a mathematical model to study LV 
remodeling post-MI by quantifying the dynamic balance between ECM construction and destruction. The 
mathematical model incorporated the key factors and demonstrated that TGF-pi stimuli and MMP-9 interventions 
with different strengths and intervention times lead to different LV remodeling outcomes. The predictions of the 
mathematical model fell within the range of experimental measurements for these interventions, providing 
validation for the model. 

Conclusions: In conclusion, our results demonstrated that the balance between ECM synthesis and degradation, 
controlled by interactions of specific key factors, determines the LV remodeling outcomes. Our mathematical 
model, based on the balance between ECM construction and destruction, provides a useful tool for studying the 
regulatory mechanisms and for predicting LV remodeling outcomes. 



Background 

Myocardial infarction (MI) is a leading cause of conges- 
tive heart failure (CHF) [1,2]. In response to the MI sti- 
mulus, the left ventricle (LV) undergoes structural and 
functional adaptations that collectively have been termed 
as LV remodeling [3]. Adverse LV remodeling pro- 
gresses to CHF in about 25% of post-MI patients, but 
the mechanisms that drive this progression remain 
poorly understood. During LV remodeling, both extra- 
cellular matrix (ECM) degradation and synthesis 
increase [4]. When ECM degradation dominates over 
synthesis, LV rupture can occur. When ECM synthesis 
dominates over degradation rates, fibrosis can occur. 
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Fibrosis increases myocardial stiffness and further 
depresses LV function to culminate in CHF [5,6]. There- 
fore, understanding what regulates the balance between 
ECM degradation and synthesis post-MI is critical to 
understand the mechanisms of LV remodeling and may 
allow us to target specific early diagnostic indicators to 
better guide treatment protocols. 

Previous studies have shown that matrix metallopro- 
teinases (MMPs) regulate ECM degradation and fibro- 
blasts regulate ECM synthesis [7-9]. MMP-9, 
transforming growth factor-Pi (TGF-Pi), tissue inhibitor 
of metalloproteinase-1 (TIMP-1), and collagen I levels 
are significantly elevated from day 1 to day 7 post-MI 
[4,10,11]. These increases are concomitant with 
increased infiltration of macrophages and activation of 
fibroblasts [12]. LV remodeling is a complex process 
that involves the spatiotemporal interactions among 
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many biological components that remains poorly under- 
stood, in part due to the lack of complete sets of experi- 
mental data and computational models. Therefore, the 
objectives of this study were to 1) identify candidate bio- 
markers of LV remodeling post-MI from ECM gene 
expression and plasma analyte analyses, and 2) establish 
a mathematical model that incorporates experimental 
results to predict LV remodeling outcomes following 
different interventions. This model would provide a tool 
to elucidate LV regulatory mechanisms, estimate un- 
measurable variables, and predict outcomes following 
multiple therapeutic scenarios. 

Results 

Identifying Key Factors 

The key factors were pre-targeted by examining the 
most significant changes in ECM gene expression in the 
infarct region at day 7 post-MI, compared to gene 
expression in the remote non-infarcted region of the 
same LV and in the LV from control group. In the 
ECM gene array analysis, total RNA yield was 1.0 ± 0.1, 
1.9 ± 0.2, and 2.9 ± 0.3 (ig/mg LV tissue for control, 
remote, and infarct samples, respectively (p < 0.05 for 
control vs remote and infarct, and for remote vs infarct). 
Of the 84 genes examined, 51 genes were differentially 
expressed among control, remote, and infarcted groups 
(all p < 0.05). The most prevalent pattern of gene 
expression changes was an increased expression level in 
the infarct tissue, compared to both control and remote 
groups. Of the 51 genes, 17 genes showed > 2.5-fold 
change in the infarct region, and these genes are listed 
in Table 1. Of the 17 genes with >2. 5-fold change, the 
most significantly over expressed genes are cadherin 3, 
collagen 1, collagen 2, and collagen 3, osteopontin, peri- 
ostin, tissue inhibitor of metalloproteinase-1, fibronectin, 
secreted protein acidic and rich in cysteine (SPARC), 
and transforming growth factor-p. From this list, col- 
lagen I was selected as a key factor, because it is the 
major collagen in the normal LV, accounting for 90% of 
cardiac ECM. Osteopontin expression levels increased 
206-fold in the infarct region compared to the control 
group, suggesting strong macrophage activation at day 7 
post-MI. TIMP-1 expressions increased 31-fold in the 
infarct compared to the control group, suggesting strong 
inhibition of proteolytic activity. Periostin gene expres- 
sion level was increased 5.5-fold, and TGF-pi gene 
expression level was increased 2.6-fold in the infarct 
sample at day 7, suggesting significant fibroblast func- 
tions post-MI [13-16]. 

We further examined plasma profiles at day 7 post-MI 
to determine if gene level changes were mirrored in the 
post-MI plasma. In the plasma profiles, 11 proteins 
increased in the day 7 post-MI samples (Table 2). Of 
these, increased MMP-9 and TIMP-1 indicate changes 



Table 1 ECM Gene Array (Data are Mean ± SD normalized 
levels.) 



Gene 


Control 




Remote 




Infarct 


Cdh3 


1 .OE-05 ± 


9.8E-06 


7.4E-05 ± 


5.5E-05 


2.8E-04 ± 1.3E-04 


Collal 


1.5E-02 ± 


2.6E-03 


1.9E-01 ± 


l.OE-01 


8.2E-01 ± 3.4E-01 


Col2al 


3.1E-06 ± 


6.3E-07 


l.lE-05 ± 


l.lE-05 


2.3E-04 ± 1.8E-04 


CoBal 


5.0E-03 ± 


1.5E-03 


8.0E-02 ± 


5.5E-02 


2.3E-01 ± 9.7E-02 


ColSal 


8.0E-3 ± 2.0E-3 


2.3E-2 ± ■ 


1.1 E-2 


4.1 E-2 ± 1.6E-2 


Ctgf 


2.5E-02 ± 


4.3E-03 


7.6E-02 ± 


5.0E-02 


1.7E-01 ± 9.1E-02 


Fnl 


9.0E-03 ± 


2.5E-03 


4.3E-02 ± 


2.7E-02 


1.6E-01 ± 6.9E-02 


Mmp2 


3.5E-03 ± 


1.2E-03 


1 .2E-02 ± 


7.9E-03 


1.8E-02 ± 7.5E-03 


Mmpl4 


1 .9E-03 ± 


4.6E-04 


6.8E-03 ± 


3.2E-03 


1.5E-02 ± 6.8E-03 


Ncaml 


6.2E-04 ± 


1.3E-04 


1 .8E-03 ± 


1.3E-03 


4.1E-03 ± 1.7E-03 


Postn 


5.3E-03 ± 


1.6E-03 


1.5E-01 ± 


9.2E-02 


2.9E-01 ± 1.2E-01 


Sparc 


1.3E-02 ± 


3.7E-03 


7.0E-02 ± 


3.7E-02 


l.lE-01 ± 4.5E-02 


Sppl 


1.3E-05 ± 


5.7E-06 


3.1E-04 ± 


4.5E-04 


2.6E-03 ± 2.2E-03 


Tgfbi 


2.0E-03 ± 


3.3E-04 


3.3E-03 ± 


1 .2E-03 


7.1E-03 ± 3.6E-03 


Thbsl 


7.2E-03 ± 


6.7E-03 


1 .2E-02 ± 


7.6E-03 


4.0E-02 ± 1 .6E-02 


Thbs2 


9.0E-03 ± 


2.3E-03 


l.lE-02 ± 


3.6E-03 


3.9E-02 ± 1 .4E-02 


Timpl 


2.0E-04 ± 


9.4E-05 


4.1E-03 ± 


3.1E-03 


6.1E-03 ± 3.3E-03 



Control is unoperated mice; Ml is 7d post-MI mice. 



in ECM remodeling. There were several proteins that 
influence macrophage activation and function, including 
macrophage inflammatory protein-1 and osteopontin 
that increased in post-MI plasma. Via analysis of the 
plasma changes post-MI [4,7-12], we further confirmed 
key factors involved in macrophage and fibroblast func- 
tions, namely collagen I, MMP-9, and TIMP-1. 

In addition, we examined the ECM productions in iso- 
lated cardiac fibroblasts stimulated with TGE-P^. Fibro- 
blast ECM array analysis showed that TGF-Pi 
stimulation of cardiac fibroblasts up regulated 5 genes 
and down regulated 7 genes, which are shown in Table 
3. Collagen I expression and SPARC expression levels 
were doubled, and TIMP-1 expression increased 5-fold 
in response to TGF-Pi stimulation at the concentration 
of 10 ng/mL. These experimental data indicated a pri- 
mary regulatory effect of TGF-Pi on fibroblast ECM 
production. Interestingly, the 5 up regulated genes are 
among the 17 ECM genes that were significantly 
expressed in the LV, indicating that the cardiac fibro- 
blast is likely the major tissue source for these ECM 
genes (Table 1). 

We also examined the correlations between LV wall 
thickness with 6 genes that were over expressed post- 
MI. The values were 0.76 for collagen lal, 0.64 for 
collagen 2a 1, 0.75 for collagen 5a 1, 0.60 for periostin, 
0.61 for osteopontin, and 0.63 for TGF-Pi. 

In summary, the key factors identified were macro- 
phages, fibroblasts, TGF-Pi, MMP-9, and collagen. 
Based on these experimental results, we developed a 
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Table 2 Multi-analyte Profiling of Control and 7 day Post-MI Plasma 





Control n = 6 


7 d Ml n = 7 


P value 


Clusterin (|jg/mL) 


330 ± 40 


510 : 


± 200 


0.046 


Cystatin-C (ng/mL) 


360 ± 20 


530 : 


± 100 


0.008 


Eotaxin (pg/mL) 


1 rr\r\ i nr\r\ 

1 500 ± 200 


1 QOO 


-1- /inn 


0.049 


riUllllULjcll \,lllLj/lllLj 


19 + 9 
1 Z X z 


19 


± 5 


n ni 9 

U.U 1 z 


Haptoglobin (|jg/mL) 


83 ± 30 


180 


± 50 


0.002 


Macrophage Inflammatory Protein-1 gamma (ng/mL) 


24 ± 4 


33 


± 9 


0.042 


Matrix Metalloproteinase-9 (ng/mL) 


71 ± 9 


96 : 


± 20 


0.009 


Myeloperoxidase (ng/(mL) 


58 ± 10 


100 


± 30 


0.006 


Osteopontin (ng/mL) 


250 ± 60 


390 : 


± 100 


0.041 


Serum Amyloid Protein (|jg/mL) 


21 ± 3 


33 


± 9 


0.013 


TIMP-1 (ng/mL) 


0.8 ± 0.1 


3.7 : 


± 2.0 


0.002 



framework of the interaction loops among the identified 
key factors (Figure 1). 

Unking Experimental Results to Mathematical Modeling 
Framework 

In the ECM construction pathway, collagen is secreted 
by fibroblasts. Growth and secretion of fibroblasts are 
stimulated by TGF-Pi (Table 3). In ECM destruction 
pathway, MMP-9 is a key factor by cleaving collagen. 
The major source of MMP-9 is the macrophage, which 
infiltrates into the infarct region post-MI. The majority 
of macrophages in the LV post-MI are differentiated 
from peripheral blood monocytes stimulated by che- 
moattractants including TGF-Pi [17]. The major source 
of TGF'Pi is the macrophage. Meanwhile, there are 



Table 3 Fibroblast ECM array In serum free control and 
10 ng/ml TGF-p stimulated fibroblasts 





Serum Free 


TGF-P 


p Value 


ECM/Growth Factors 


Collal 


4.842 ± 1 .399 


9.614 ± 3.324 


0.028 


Col5al 


0.301 ± 0.085 


0.605 ± 0.230 


0.038 


Fblnl 


0.012 ± 0.005 


0.008 ± 0.004 


0.032 


Sparc 


2.709 ± 0.204 


4.932 ± 0.379 


<0.001 


Tgfbi 


0.022 ± 0.016 


0.013 ± 0.015 


0.009 


Cell Adhesion Molecules 


Ncaml 


0.043 ± 0.016 


0.169 ± 0.054 


0.011 


Pecaml 


0.000060 ± 0.000012 


0.000022 ± 0.00001 1 


0.015 


Sgce 


0.187 ± 0.035 


0.146 ± 0.018 


0.030 


Vcaml 


0.301 ± 0.136 


0.161 ± 0.088 


0.031 


MMPs^lMPs 


Mmp7 


0.0000075 ± 0.000001 1 


0.0000056 ± 0.0000004 


0.026 


Timpl 


0.052 ± 0.022 


0.272 ± 0.136 


0.035 


Timp2 


0.326 ± 0.109 


0.188 ± 0.050 


0.027 



Data are Avg ± SEM levels (2"^'"^) normalized to 5 housekeeping genes for n = 
4 paired fibroblast sets. 



interactions between the ECM construction and destruc- 
tion pathways: a) MMP-9 regulates ECM construction 
by activating TGF-P i which stimulates collagen synth- 
esis; b) TGF-P 1 induces TIMP-1, which inhibits collagen 
degradation by blocking MMP-9 activity. Linking these 
key factors with their sources and effects allows us to 
develop the mathematical model for quantitative 
analysis. 

Mathematical Modeling 

We established a set of nonlinear differential equations 
to model the temporal interactions among the key fac- 
tors identified by our experimental results. The model 
incorporated the following variables: macrophage cell 
density (Mo, cells/mm^), fibroblast cell density (F, cells/ 
mm^), collagen concentration (C, (ig/(iL), TGF-P i con- 
centration (Tp, pg/(iL), and activated MMP-9 concentra- 
tion (M9A, pg/(iL). Rates of cell number change were 
determined by the summation of constructive effects 
(migration rate or proliferation rate) and destructive 
effects (death rate or removal rate). Rates of chemical 
factors (TGF-P 1, MMP-9, collagen, etc) change were 
determined by the net difference between the synthesis 
rate and degradation rate. 

Four assumptions were used based on experimental 
results: 1) All monocytes that migrate to the infarct 
region are differentiated to macrophages [18]; all acti- 
vated macrophages are differentiated from peripheral 
blood monocytes since previous studies have shown that 
<5% of macrophages undergo mitotic division [19]; 2) 
The major source of fibroblasts is the resident cell and 
the contribution of circulating fibroblasts is ignored 
[20]; 3) The majority (80%) of TGF-Pi secreted at the 
injured site becomes activated [21,22]; 4) ECM proteins 
in the infarct region are secreted by local cells. Our 
experimental data showed that ECM gene expression 
was higher in the infarct region than the expression in 
the remote region, and ECM gene levels in the remote 
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Figure 1 A working model of post-MI scar formation and remodeling with constructive and destructive modules. Red arrows represent 
the destruction pathway and green arrows represent the construction pathway. The regulation scheme includes lines in black, which can be a 
driving or a triggering stimulus. LV remodeling outcomes are denoted as collagen concentrations in the scar tissue. 

V J 



region were higher than that in the LV of controls, indi- 
cating a local source of ECM production. 

Accordingly, the scar formation post-MI was modeled 
by the following set of differential equations 



inf iltration removal 

^crowding = kmemit) + J^^/PM^ +F/pf + C/ pc, 

fihrohlast proliferation rate with crowding effect apoptosis 
F' = kpF^iTf^) F[l - kcrowding] - dpF , 

fibroblast tcfB 
secretion secretion activation ^^S^adation 



M9A + C 



koff 



[CM9 



CID + MMP9 



(1) 

(2) 
(3) 

(4) 
(5a) 



secretion by activation by TGF-p ocupiedfor collagen degradation inactivation 



h[T) 



- koff[CMs] - kon^SAC -dM9M^A, 



collagen synthesisby fibroblast sitmulated by TGF—fiwith collagen crowding collagen degradation by M M P9 



+ koff [CMo] - kon M<)A C, 



(6) 



(7) 



The parameters used in these equations with their 
biological meanings, experimental values, units, and 
references were listed in Table 4. The interaction func- 
tions M(Tp), Fg{Tp), and Fc(Tp), were established based 
on in vivo, and in vitro experiments [17,23-26]. We 
employed the function Fc(C) in a form given by Waugh 
and colleagues [27,28]. The forms of these functions we 



took in this study were described in equation 8. 

M{Tfi) = 0.3350T^ - 6.3095T^ + 32.2810T^ + 57.3019 
Fg(T^) = 0.0492T^ - 0.9868T^ + 6.5408T^ + 7.1092 
Fc{Tfi) = 0.0092T^ - 0.1552T^ + 0.6279T^ + 0.2527 
Fc{C) = -4.33 X 10-^°C^ + 0.0000009C^ -0.00055C + 0.13 



(8) 



Plot of these constructed functions and the available 
experimental data were shown in Figure 2. 

Equation 1 determines the rate of macrophage (Mo) 
infiltration. The function M(Tp) describes the migration 
rate of macrophages to the scar tissue [17]. Since we 
assumed that all monocytes differentiated into macro- 
phages (assumption 1), M(Tp) also represents the migra- 
tion of monocytes stimulated by TGF-Pi. Parameter 
djvio denotes the emigration rate of macrophages [29], 
since macrophages do not die locally in the scar tissue 
but emigrate to the lymph node system for disposal. 

Equation 2 determines the crowding effect of myo- 
cytes, endothelial cells, vascular smooth muscles cells, 
fibroblasts, macrophages, and collagen in the myocar- 
dium, which are affected by total environment density 
(5b) [30,31]. Parameter kmemit) = 0.17 + 0.56^"^-^^ denotes 
the crowding effects contributed by myocytes, endothe- 
lial cells, and vascular smooth muscle cells. The para- 
meter 0.17 represents the total percentage of endothelial 
and smooth muscle cells that account for 7% and 10% 
of total cell numbers in normal mouse myocardium, 
respectively [32]. Parameter 0.56 represents the percen- 
tage of myocytes (56%) in normal mouse myocardium 
[32]. Parameter dmc = 0.05 represents the rate of myo- 
cyte cell death since ischemic myocytes undergo necro- 
sis in the infarct region post-MI. The crowding effect of 
macrophages, fibroblasts, and collagen was considered 
by calculating the normalized density with respect to 
their maximum density in scar tissue [28]. Temporal 
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Table 4 Pre-determined parameters in the mathematical model 



Symbol 


Biological meaning 


Value 


Units 


Ret 




Macrophage removal rate**(eqn 1) 


0.6 


day 


[29] 


PMO) 


maximal macrophage density (eqn 2) 


2500 


cells/mm^ 


[28] 


PF 


maximal fibroblast density(eqn 2) 


1 250 


cells/mm^ 


[28] 


Pc 


maximal collagen density (eqn 2) 


3300 


|jg/mm^ 


[28] 


kp 


Fibroblast growth rate* (eqn 3) 


0.924 


day"^ 


[62] 


dp 


Fibroblast apoptosis rate (eqn 3) 


0.1 2 


day"^ 


[33] 




Macrophage TGF-p production rate (eqn 4) 


0.07 


pg/cell/day 


[34] 




riUlUUIdbL 1 Vjr p piUUULLIUII IdLc \,cqil ^) 




py/Lcii/udy 


pjj 


dip 


TGF-(3 degradation rate^(eqn 4) 


15 


day"^ 


[36] 




Macrophage secretion MMP9 rate (eqn 5) 


3 


pg/cell/day 


estimated 




MMP-9 degradation rate(eqn 5) 


0.875 


day"^ 


[63-65] 




Kinetic reaction speed (eqn 5) 


3 X 10"^ 


l/(|jg/mm^)s"^ 


[42,66,67] 


koff 


Kinetic reaction speed (eqn 5) 


4 X 10"^ 


s-^ 


[42,66,67] 


konc 


Kinetic reaction speed (eqn 6) 


0.004 


s-^ 


[42,66,67] 


kpc 


Fibroblast collagen production rate (eqn 6) 


20 


|jg/cell/day 


[41] 



*The growth rate of cells was determined by the population doubling time (T2) via equation k = In 2/T2 [68]. 



Since macrophages emigrate from the scar tissue to lymph node system instead of dying locally in the scar tissue, the removal rate of macrophage, dMA/ was 
used in our model. 

¥ The decaying rate of chemical factors was calculated from their half-life (T1/2) via the equation d = In 2/T1/2 [68]. 



profiles of total crowding effects, crowding effects con- 
tributed by macrophages, fibroblasts, and collagen, and 
the parameter kymm (^) were shown in Figure 3. 

Equation 3 determines the rate of fibroblast (F) den- 
sity changes based on the assumption that majority of 
fibroblasts come from the proliferation of resident cells 
(assumption 2). The function Fg{Tp) denotes the stimu- 
lating effects of TGF-Pi on the growth rate of fibroblasts 
[17]. Parameter dp represents the apoptosis rate of the 
fibroblast [33]. 

Equation 4 determines the rate of TGF-Pi concentra- 
tion change, wherein kpx denotes the TGF-Pi secretion 
rate of fibroblasts [34] and kjvioT denotes the TGF-Pi 
secretion rate of macrophages [35] since a major source 
of TGF-Pi in the scar tissue is the activated macro- 
phage. Parameter dxp represents the degradation rate of 
TGF-Pi, which can be calculated from the half life data 
[36]. 

TGF-Pi gene levels demonstrated temporal progres- 
sion at the early stage post-MI. Gene expression profile 
of TGF-Pi increased post-MI, peaked at day 2, and 
returned to normal levels after day 7 in mice post-MI 
[37]. Since majority of TGF-Pi secreted in the infarct is 
activated (Assumption 3), gene expression profile can be 
used as an activation pattern of TGF-Pi. The function, 
ut denotes temporal profile of TGF-Pi activation post- 
MI [38-40] and works as the input in the simulation. It's 
worth to mention that TGF-pi activation levels were 
different with respect to different scar size. A temporal 



profile of ut representing small scar was shown in Fig- 
ure 3. 

Equation 5 determines the rate of activated MMP-9 
concentration change. Proteolytic collagen degradation 
with activated MMP-9 is described in equation 5a, 
where M9A, C, CM9, and CID denote activated MMP-9, 
collagen, binding of MMP-9 and collagen, and degraded 
collagen peptide concentration, respectively. MMP-9 is 
inhibited primarily by TIMP-1, and TIMP-1 is induced 
by TGF-Pi. Thus, we established an inhibition function 
h(T) = l/(l+Tp/TpN) with TpN = 6.0 pg/(iL to represent 
the inhibition effect. 

Equation 6 determines the rate of collagen concentra- 
tion changes. Collagen secretion rate by fibroblasts was 
denoted by parameter/C^c- Meanwhile, the function, 
(Tp), characterizes effects of TGF-pi on collagen secre- 
tion rate by fibroblasts [23,41]. Function, Fdc), denotes 
the effect of collagen density on fibroblast secretion rate 
[27,28]. 

Equation 7 determines the concentration change of 
CM9, based on the theoretical model for collagen degra- 
dation by MMPs proposed by PopeFs group [42,43]. 

Computational simulations 

Computational simulations of scar formation (collagen 
deposition) were carried out by solving the nonlinear 
differential equations with MATLAB. Initial conditions 
of the fibroblast and macrophage densities were chosen 
as F(0) = 20 cells/mm^ and Mo(0) = 5 cells/mm^ 
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150 




0 2 4 6 8 10 
TGFp concentration (pg/mL) 



D) 




0 2 4 6 8 10 
TGFp concentration (ng/nriL) 




0.12 



O 




0 2 4 6 8 10 0 500 1000 1500 

TGFp concentration (ng/nnL) Collagen density (jli g/mL) 

Figure 2 Macrophage migration rate IVI(Tp) [17], fibroblast growth rate Fg(Tp) [17], and fibroblast secretion rate Fc(Tp) [23,41] plotted as 
functions of TGF-Pi concentration. The original experimental data (x. Mean ± SD) are shown for comparison. The function Fc(C) was provided 
by Waugh et al [27,28]. 



Accordingly, Tp(0) = 0.21 pg/(iL, M9a(0) = 7.1 pg/(iL, C 
(0) = 839.5 [ig/[iU CM9(0) = 447.6 [ig/[iL were calcu- 
lated by the equilibriums of equations 1-7 for normal 
LV. All the simulations shown in this study used the 
same initial conditions. The initial conditions were cho- 
sen based on measurements in the normal LV for both 
the control and MI groups (MI induced at day 0). The 
simulations covered the LV remodeling process from 
day 0 to day 30 post-MI. 

Model validation 

To validate our mathematical model, we compared our 
simulation results to experimental data from our lab or 
reported in the literature [9,12,44]. We normalized 
macrophage and fibroblast cell densities and MMP-9 
concentrations to the corresponding measurements in 
the normal LV and plot the experimental data (Mean ± 



SD) in Figure 4. Our computational simulation results 
were also normalized to the corresponding initial condi- 
tions. This normalization will give us the relative fold 
changes of macrophage and fibroblast cell densities and 
MMP-9 concentration. The simulation results showed 
similar progression trend of cell densities and MMP-9 
concentration profiles, peak values, and stable values as 
the experimental results (Figure 4). 

In addition, our simulations correctly predicted MMP- 
9 responses to three TGF-Pi stimuli corresponding to 
reduced, normal, and elevated post-MI activation 
strength. Others have reported an early increase of 
MMP-9 levels of 78 ± 19 pg/(iL for small infarcts and 
195 ± 63 pg/(iL for large infarcts [45,46]. These experi- 
mental results agreed with the MMP-9 predictions by 
our mathematical model, which peak at 80 pg/(iL for 
small infarcts, and 220 pg/(iL for large infarcts. 
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10 20 
Time 



10 20 
Time 



30 





05 

o 0.95 
I- 

0.94 

0 10 20 30 0 10 20 

Time Time 

Figure 3 Plots of crowding effects (equation 2) and TGF-pi activation function Uj (equation 4) A: Crowding effects [kmem (t) in equation 
2) contributed by myocytes, endotlielial cells, and vascular smooth muscle cells. B: Crowding effects contributed by macrophages, fibroblasts, 
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respectively. MMP-9 profiles were shown in Figure 5 
and details of the simulation setup for Figure 5 were 
explained as follows. 
Effects of TGF-pj levels 

We employed activation of TGF-Pi, Ut in equation 4, at 
reduced, normal, and elevated post-MI strength. The 
activation function peaked at 15, 30, and 60 pg/(iL/day, 
according to active expression levels observed in small, 
median, and very large infarcts, respectively. Temporal 
profiles of cell densities of macrophages, fibroblasts, 
concentrations of MMP-9, collagen, and TGF-Pi 
responses to the stimuU Ut are shown in Figure 5. For 
the TGF'Pi stimulus level setting at 30 pg/(iL/day (a 
normal post-MI level in mice), TGF-Pi peaked after day 
2 post-MI, macrophage density peaked at day 3, MMP-9 
concentration peaked at day 4, and all were returned to 



normal levels at 30 days post-MI. In contrast, fibroblast 
density and collagen density continued to increase 
beginning at day 4, reached a stable value after day 20 
post-MI, and then remained at a higher equilibrium 
level at day 30 post-MI (blue solid line in Figure 5). 
This computational predictions agreed with the experi- 
mental observations on the peak time of macrophages 
[9,12], and progression trend of MMP-9 [10] and col- 
lagen for stable LV remodeling [47]. 

In the case of reduced TGF-Pi levels (with an activa- 
tion peak at 15 pg/(iL/day in Figure 3), less monocytes 
were attracted to the infarct. Thereby, less macrophages 
and MMP-9 appeared at the early stage. The reduced 
MMP-9 levels slowed down collagen degradation, lead- 
ing to a higher collagen concentration (black dotted 
lines in Figure 5. This simulation agreed with the 
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temporal profiles of macrophage and fibroblast density 
[9,12], and relative low expression levels of MMP-9 for 
moderate size of infarcts compared to large infarcts 
[46,48]. 

When the TGF-Pi stimulus strength was elevated to a 
level 2-fold higher than normally seen post-MI in mice, 
more macrophages infiltrated to the infracted region in 
the early days (Figure 5 red dash lines). Elevated macro- 
phage infiltration led to high levels of crowding effect 
and high concentrations of MMP-9 for collagen degra- 
dation. Therefore, fibroblast growth was inhibited earlier 
due to the crowding effect, which led to less collagen 
secretion and a net negative collagen deposition. These 
results indicated an increased early susceptibility to LV 
rupture between day 2 to 6 post-MI and prolonged LV 
remodeling for those mice that did not rupture. Though 
we know that LV ruptures frequently occurred at day 2- 
4 post-MI in a normal remodeling process [49], this 
simulation has not been examined experimentally and 
raises some interesting hypotheses. 
Effects of MMP-9 interventions on ECM destruction 
With the validated parameter settings, we also used the 
mathematical model to predict the effects of MMP-9 
interventions at different strengths and intervention 
times. Specifically, we simulated the LV remodeling 
responses to three different MMP-9 interventions post- 
MI: 1) elevation of MMP-9 level (200 pg/(iL) beginning 
at 8 hours post-MI to mimic the earlier increase of 
MMP-9 levels seen with reperfusion, 2) elevation of 
MMP-9 levels (200 pg/(iL) beginning at 7 days post-MI 
to mimic a prolonged macrophage infiltration, and 3) 
reduced elevation of MMP-9 levels (100 pg/(iL) begin- 
ning at 7 day to mimic therapeutic targeting of MMP-9. 
The LV remodeling responses were shown in Figure 6. 
Early MMP-9 intervention (scenario 1) significantly 
decreased collagen density at day 30 post-MI, consistent 
with the beneficial effects seen with reperfusion. Ele- 
vated MMP-9 levels (scenarios 1, 2, and 3) led to 
reduced collagen density post-MI, compared to collagen 
levels seen in LV remodeling outcomes without MMP-9 
intervention, when MMP-9 levels would be starting to 
fall by day 7. Compared to scenario 2, reducing MMP-9 
intervention levels at day 7 post-MI (scenario 3) had 
less effect on collagen degradation, suggesting that the 
primary effects of MMP-9 might occur before day 7. 

Discussion 

This study is the first investigation to integrate mathe- 
matical modeling with ECM and fibroblast gene array 
data and plasma analytes to predict ECM remodeling 
post-MI. We have integrated in vivo, in vitro, and in 
silico approaches to dissect the complicated interactions 
among multiple regulatory factors in LV remodeling. As 
such, this study provides a promising modeling 



approach for elucidating the complicated LV remodeling 
process. The most important findings of this study are 
highlighted as follows. 1) The balance between ECM 
construction and destruction kinetics is likely the key 
determinant of scar formation, and interventions to 
either the construction or destruction side lead to differ- 
ent remodeling outcomes. 2) Dynamic interactions 
among key factors in LV remodeling determine LV 
remodeling outcomes post-MI. Altering initial condi- 
tions, intervention strengths, or intervention times have 
significant effects on LV remodeling outcomes, and 
these effects could be simulated accurately by our 
model. 3) Collagen, MMP-9, TGF-Pi, and TIMP-1 are 
critical biomarker candidates of LV remodeling 
outcomes. 

Our experimental results on microarray and plasma 
data provided the foundations to build our computa- 
tional framework. We examined 84 ECM genes and 
chose the 17 genes that were most highly over- 
expressed in the infarct region compared to both con- 
trol and non-infarcted groups (>2. 5-fold over-expres- 
sion). The expression levels of several of these factors 
were further verified by our plasma data at protein level. 
One interesting finding was that SPPl (osteopontin) 
gene expression levels increased 206-fold in the infarct 
region compared to the control group at day 7 post-MI, 
suggesting strong macrophage activation. The plasma 
protein levels of osteopontin increased from 250 ± 60 
ng/mL in controls to 390 ± 100 ng/mL in post-MI sam- 
ples, which adds support for the critical role of macro- 
phages in our mathematical model. Therefore, the 
primary selection of the most highly changed genes 
allowed us to focus on the most significant factors at 
gene level and predict the possible interactions at pro- 
tein level and cellular level. 

An interesting observation was that MMP-9 mRNA 
levels did not increase in our gene array analysis, but 
MMP-9 protein levels increased in the plasma data. It is 
well known that MMP-9 protein levels and activation 
were increased at day 7 post-MI. This conflicting phe- 
nomena is caused by pre-formed MMP-9 proteins 
stored in leukocytes, which do not rely on increased 
gene expression to induce MMP-9 secretion and activa- 
tion [50]. In addition, the shift from a normal ventricle 
composed of cardiac myocytes, endothelial cells, and 
smooth muscle cells to an infarcted ventricle comprised 
primarily of cardiac fibroblasts and inflammatory cells 
could result in a quantitative no gain in gene expression 
level but qualitative increase in MMP-9 function due to 
increased activation and increased substrate availability. 
Karl Weber's group showed that MMP-9 protein levels 
increased early post-MI, but the mRNA levels were not 
elevated [51]. We have also shown that MMP-9 activa- 
tion increases upon reperfusion in a dog model of 
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cardiac ischemia/reperfusion [7]. These results support 
our observations in this study and highlight the fact that 
MMP-9 is primarily regulated at the post-translational 
level Furthermore, the fact that MMP-9 gene levels do 
not change at day 7 post-MI also suggest that MMP-9 
should be evaluated at the protein level 

There is limited data available to construct suitable 
functions of M{Tp), Fg{Tp), F,{Tp), and F,{c), Wahl et al 
presented a set of in vitro data on monocyte chemotaxis 
induced by TGF-P with concentration ranges from 0.05 
pg/mL to 10 pg/mL [17]. Their experiments also pro- 
vided production of fibroblast growth activity induced 
with TGF-p concentration range from 10 pg/mL to 10 
ng/mL. These data were used in this study to construct 
the functions of M {Tp) and Fg{Tp), Roberts et al have 
presented a set of data on in vitro collagen formation 
stimulated by TGF-P with the concentration range from 
0.0825 pg/mL to 2.5 pg/mL [23]. Ignotz and Massague 
have examined the effects of TGF-p on coUagenases in 
chick embryo fibroblasts stimulated by TGF-P at a con- 
centration range of 1.25 ng/mL to 25 ng/mL [41]. We 
have done experiments to quantify the ECM production 
of cardiac fibroblast stimulated with TGF-p (10 ng/mL). 
These experimental data were used to build and justify 
the function, FdTp), Loftis and colleagues have studied 



effects of collagen density on cardiac fibroblast behavior 
and showed elevated fibroblast activities stimulated with 
higher collagen concentrations (750 (ig/mL - 1250 \ig/ 
mL) [52]. Wu and colleagues showed that propeptides 
at lower levels caused about 80% decrease in collagen 
synthesis compared to control [53]. In this study, we 
used an in silico function Fdc) provided by Waugh and 
colleagues [27,28]. These experiments were done by dif- 
ferent groups with different conditions and the in vitro 
experiments might not reflect the in vivo interactions. 
In addition, we established these functions using polyno- 
mial interpolation. It's likely that these functions are not 
the optimal forms to describe the in vivo interactions; 
however, they provide us a baseline for further study. 

Our computational results demonstrated that altering 
the strength of TGF-P i altered LV remodeling out- 
comes. Elevated TGF-Pi levels at the early stage (day 3 
post-MI) led to elevated macrophage density and MMP- 
9 levels, decreased fibroblast secretion of collagen and 
collagen deposition, and thereby, prolonged the progres- 
sion of remodeling. Wetzler and colleagues have shown 
that the prolonged persistence of macrophages at the 
late phase (after day 7 post-MI) impairs the wound heal- 
ing process [54]. It has also been shown that elevated 
TGF-P 1 levels delays wound healing post-MI [55,56]. 
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Our simulation results are in agreement with these pre- 
viously published studies. 

Simulations of different TGF-Pi strengths also shed 
light on the regulation scheme of ECM construction 
and destruction. For ECM construction regulation, 
active TGF-Pi stimulated fibroblast proliferation and 
collagen secretion, which increased the crowding 
effect. The increased crowding coefficient ramped 
down fibroblast proliferation and TGF-Pi secretion 
(negative feedback), which slowed down the stimulus 
for monocytes to migrate into the infarct region. The 
decrease of monocytes number led to less macrophage 
infiltration, which then reduced the crowding coeffi- 
cient. Meanwhile, reduced macrophages lead to less 
TGF-Pi secretion by macrophage, which further slo- 
wed down collagen synthesis. For collagen destruction 
regulation, TGF-Pi induced macrophage infiltration, 
which lead to elevated MMP-9 secretion, elevated col- 
lagen degradation, and thereby reduced crowding coef- 
ficients. Smaller crowding coefficients lead to elevated 
fibroblast proliferation and collagen secretion. Notably, 
there are two types of negative feedback schemes in 
the mathematical model: degradation (apoptosis or 
emigration) rates associated with proteins (cells) and 
the crowding effects. Degradation rates are constants 
and determine how fast the proteins (cells) can 
respond to stimuli. Crowding effects are time varying 
impacts imposed by the environment. Through these 
regulation schemes, a dynamic balance of collagen 
construction and destruction can be maintained to 
generate a stable scar. Furthermore, profiles of crowd- 
ing effects elucidated the transition from the normal 
LV to scar tissue with respect to cell types and col- 
lagen concentrations (Figure 3). 

It is worth mentioning that there exist biological nega- 
tive feedbacks in our mathematical model. For macro- 
phage density regulation, there was a positive feedback 
loop containing macrophage and TCF-P^: TGF-Pi sti- 
mulated monocytes migration, leading to macrophage 
infiltration; macrophages secreted TGF-Pi which might 
attract more macrophages to the infarct site. We 
observed elevated macrophage density and MMP-9 con- 
centrations corresponding to increased TGF-Pi levels 
(Figure 5). This positive feedback loop was inhibited by 
emigration of macrophages and degradation of TGF-Pi. 
In addition, Wahl et al pointed out that the monocyte 
chemotactic activity increased in response to low con- 
centrations of TGF-Pi stimuli, while the chemotactic 
activity decreased in response to higher concentrations 
of TGF-Pi [17]. Therefore, as TGF-Pi levels continu- 
ously increased, infiltration speed of macrophages 
decreased as shown in figure 2 (M(Tp)), suggesting a 
secondary biological inhibition scheme of the TGF-Pi - 
macrophage positive feedback loop. 



There are a few limitations of the mathematical model 
that resulted to a large degree from the model assump- 
tions. More research is needed to address these Umita- 
tions and further enhance the models. First, our model 
calls for accurate determination of MMP-9 activation 
and inhibition functions. Complete time-course mea- 
surements of TIMP-1 and the other three TIMPs (for 
MMP-9 inhibition) and MMP-3 (for MMP-9 activation) 
would provide additional details on the regulators of 
MMP-9 function. Second, the large differences between 
simulation results and experimental measurement of 
MMP-9 concentrations before day 4 post-MI in Figure 
4 is likely due to the effect of neutrophils, an early 
source of MMP-9, indicating the need to investigate the 
role of neutrophils in the early stage (days 1-5) post-MI 
[7]. This inconsistency of the simulated MMP-9 concen- 
trations and measurements from experimental results is 
one of the approaches for us to employ more biomar- 
kers and key factors in model development. Third, 
Interleukin IL-1, IL-6, IL-10, and tumor necrosis factor- 
a profiles need to be included to better study macro- 
phage activation phenotypes [57]. Fourth, the effects of 
other chemoattractants on monocyte migration on the 
LV remodeling need further investigation. The predicted 
fibroblast density at days 14-30 post-MI was lower than 
experimental measurements. This might be caused by 
assumption 2, that the increase of fibroblast numbers 
post-MI was contributed solely by proliferation of local 
cells. We will investigate the migration rates of fibro- 
blasts in the future to refine our mathematical model. 
Regardless of the above limitations, we were able to use 
our model to compute the progression of macrophages, 
fibroblasts, and collagen density, MMP-9 and TGF-Pi 
concentrations post-MI. 

This systems biology study for LV remodeling can be 
expanded to include proteomics and cardiac functions 
in future studies. We employed plasma data in this 
study since plasma proteins reflect the process of LV 
remodeling and plasma data are more directly translata- 
ble to the clinic. However, we are well aware that mea- 
suring tissue protein levels will provide a more direct 
evaluation of LV remodeling. Further investigation on 
the ECM proteomics in cardiac samples has been 
planned in our future research to establish a more com- 
plete mathematical model for LV remodeling. Though it 
is beyond the scope of the current paper, we have pre- 
viously reported some data on cardiac function [58-60] 
and our future study will integrate cardiac function into 
the mathematical model. 

The computational model for post-MI LV remodeUng 
developed here illustrated the dynamic interactions 
among critical factors in LV remodeling. This is the first 
mathematical model focusing on the protein and cellular 
interactions post-MI. Thus, this model provides a strong 
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foundation for future studies to build a more compre- 
hensive model that takes into account a more complete 
set of parameters. The model also provides a tool to 
guide experimental designs by identifying candidate fac- 
tors to intervene, the proper intervention time, and 
doses for effective interventions to achieve the most 
beneficial outcomes. As an example, we have shown the 
effects of MMP-9 intervention time and doses on LV 
remodeling outcomes in this study. Though this model 
was established based on in vitro data and in vivo data 
from mice, the modeling approach can be used to 
develop models for other scenarios such as the LV 
remodeling of human MI and LV remodeling under 
reperfusion conditions. 

Conclusions 

In conclusion, we developed a set of differential equa- 
tions to quantitatively model the dynamic interactions 
and temporal changes of the key components identified 
from our experimental results. Predictions of the mathe- 
matical model fell well within experimental measure- 
ments, particularly with regard to macrophage 
infiltration and matrix remodeling. This mathematical 
model provides a powerful tool to better understand 
how the dynamic balance between ECM construction 
and ECM destruction influences LV remodeling 
outcomes. 

Methods 

Mice 

All animal procedures were conducted in accordance 
with the Guide for the Care and Use of Laboratory Ani- 
mals (National Research Council, 1996) and were 
approved by the Institutional Animal Care and Use 
Committee at the University of Texas Health Science 
Center at San Antonio. Male C57BL/6J wild type adult 
mice (n = 13) at age 8.0 ± 0.5 months were used. One 
group (n = 6) served as unoperated controls, while the 
other group (n = 7) underwent coronary artery ligation 
for 7 days as described previously [9] and were sacri- 
ficed at day 7 post-myocardial infarction. 

In Vivo Procedures 

Blood was collected from the jugular vein and placed in 
a heparinized tube for plasma collection. Tissue was col- 
lected for the gene array analysis as described previously 
[9]. 

Microarray and Plasma Analysis 

Total RNA was isolated using the TRIzol plus Total 
RNA purification kit (Invitrogen). The RT^ qPCR Pri- 
mer Array for Extracellular Matrix and Adhesion Mole- 
cules (SuperArray catalog APMM-013A) was used for 
the gene array. Results were analyzed based on the 



AACt method with normalization of raw data to the 
GAPDH housekeeper gene. Data are presented as aver- 
age 2'^^^ levels. 

For the fibroblast ECM microarray analysis, cardiac 
fibroblasts were isolated from adult C57BL/6J mice and 
stimulated with or without 10 ng/ml TGF-Pi for 24 
hours [61]. Plasma samples (80 [iL) were analyzed for 67 
antigens using the quantitative immunoassay panel for 
mice provided by Rules Based Medicine in Austin, TX 
(http://www.rulesbasedmedicine.com). 

Statistical Analysis 

Control, remote, and infarct LV groups were analyzed 
by ANOVA, with the Bonferroni post hoc test. Control 
and MI plasma was analyzed by Students t-test. Unsti- 
mulated and TGF-Pi stimulated fibroblast groups were 
analyzed by paired t-test. A p < 0.05 was considered sta- 
tistically significant. 
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